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Abstract 


The  purpose  of  this  study  was  to  investigate  the  effectiveness  of  using 
fractional  order  time  derivatives  of  position  feedback  signals  in  the  active 
control  of  vibration.  To  that  end.  circuitry  was  built  and  evaluated  that  finds 
the  half  order  derivative  of  an  input  signal.  The  magnitude  response  of  the 
fractional  derivative  circuits  is  very  good,  but  there  are  large  phase  shifts 
present  that  may  degrade  the  performance  of  the  controllers.  Two  methods 
of  incorporating  the  fractional  derivative  signals  into  the  controller  were 
examined.  One  method  involved  developing  a  similarity  transformation  that 
transforms  an  integer  order  state  controller  into  an  equivalent  controller 
using  fractional  or  mixed  fractional  and  integer  order  signals.  The  second 
method  was  a  form  of  traditional  pole  placement  techniques  that  allowed  the 
direct  design  of  a  controller  using  fractional  and  integer  order  feedback 
signals.  Designs  were  tested  on  an  inverted  cantilever  beam  with  control 
provided  by  a  shaker,  and  compared  to  theoretical  predictions  and  to 
traditional  integer  order  controllers.  There  is  a  discrepancy  in  the  modeling 
that  manifests  itself  as  a  large  offset  in  magnitude,  particularly  at  very  low 
frequencies.  However,  controllers  utilizing  fractional  order  feedback 
provided  nearly  identical  control  authority  to  the  traditional  integer  order 
designs. 
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Eiperimental  Feedback  of  Fractional  Order  States 
of  a  Lightly  Damped  Structure 


I.  Introduction 


The  increasing  size  of  current  and  proposed  space  platforms  has  created 
new  problems  in  the  control  of  unwanted  vibrations.  Because  of  the  weight 
penalties  on  putting  a  platform  in  orbit,  structures  are  made  as  light-weight 
as  is  possible.  This  combination  of  increasing  size  and  minimal  weight 
construction  inevitably  leads  to  problems  with  large  amplitude,  low 
frequency  vibrations.  This  trend  runs  contrary  to  the  increasingly  stringent 
stability  requirements  needed  for  proposed  missions.  The  Hubble  Space 
Telescope  experienced  unwanted  vibrations  from  thermal  stresses  that 
interfered  with  its  mission.  The  proposed  NASA  space  station  needs  to 
remain  stable  fm*  docking  manuevers  with  the  shuttle.  The  Strategic  Defense 
Initiative  in  particular  requires  very  precise  pointing  maneuvers  and  rapid 
settling  time  after  retargeting. 

All  of  these  situations  have  led  to  an  increasing  interest  in  the  active 
control  of  vibration.  Passive  systems  provide  an  effective  means  of 
controlling  higher  frequency  vibration,  but  are  of  limited  effectiveness  in 
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controlling  the  low  frequency  modes  with  their  potentially  large 
displacements.  The  most  likely  solutions  will  involve  whatever  combination 
of  passive  and  active  damping  most  effectively  fulfill  specific  mission 
requirements.  This  thesis  will  focus  eiclusively  on  the  area  of  active  control 
of  vibration.  Extensive  research  has  taken  place  in  this  area  (See  especially 
(8)  and  (17)  which  survey  research  in  the  field),  but  many  problems  remain 
to  be  solved. 

The  current  research  effort  grew  out  of  an  interest  in  controlling 
viscoelastically  damped  structures  more  precisely.  Fractional  derivatives  of 
stress  and  strain  are  effective  in  modelling  viscoelastic  behavior  (3-7,  28), 
and  more  amenable  to  an  analytical  treatment  than  an  arbitrary  curve-fit  of 
a  damping  material's  stress-strain  curve,  klonoski  (14)  first  tested  direct 
ana)''!  feedback  of  fractional  derivatives  of  position.  Yang  (30)  developed  a 
method  of  transforming  traditional  LQR  gains  to  allow  for  the  possibility  of 
feeding  back  fractional  derivatives  of  position  and  velocity  in  a 
mathematically  optimal  fashion.  The  main  objective  of  this  thesis  is  to 
demonstrate  that  Yang's  work  on  active  damping  of  a  structure  utilizing 
fractional  derivative  feedback  is  viable  in  an  experimental  situation.  An 
alternative  approach  using  pole-placement  techniques  (13)  will  be  examined 
as  well. 

To  that  end,  circuitry  to  obtain  the  required  fractional  derivative 
feedback  was  built  and  tested,  the  experimental  structure  and  apparatus 
was  modeled  and  verified,  and  control  of  the  structure  was  attempted  for 
both  integer  and  fractional  state  feedback. 
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Fractional  calculus  is  a  little  known  cousin  of  the  classical  calculus,  and 

d''  [f(x)l 

is  nearly  as  old.  All  scientists  and  engineers  will  readily  recognize - 

dt'' 

as  the  derivative  with  respect  to  time  of  a  function  f(x),  but  very  few  have 
considered  what  would  happen  if  v  is  not  an  integer.  The  history  of  fractional 
calculus  begins  in  1693  with  L'Hospital  asking  this  question  of  Leibnitz. 
Leibnitz  said  in  his  reply  that  it  would  lead  to  a  paradox,  but  that  “Some  day 
it  would  lead  to  useful  consequences”  (23).  Euler,  Laplace,  and  Fourier  all 
considered  the  problem  at  various  times,  but  it  was  the  contributions  of 
Liouville,  Abel,  Heaviside  and  Riemann  (23)  that  finally  formalized  the 
theory  of  fractional  order  derivatives. 

A  fractional  derivative  is  a  linear  operator  that  generalizes  the  concept 
of  order  of  differentiation  from  integer  powers  only  to  fractional  values.  The 
extended  Riemann-Liouville  fractional  derivative  operator  is  defined  (2)  as: 


I  ,  ,1  dr  w(t) 
w(t)|  -  - dT  . 

‘“j  r(i-v)(t-i)'' 

0 


0  $  Vi  1 


where  T  is  the  Gamma  function  given  by 
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This  definition  is  consistent  with  the  rules  that  we  expect  of  a  derivative 
operator.  It  is  a  linear  operator  and  the  composition  rule  is  valid  (2:  8-10), 

[w(t)  +  vlt)]  -  [wit)]  +  [v(t)|  (3) 

and 

[w(t)|]  =  [w(t)l  (4) 

It  IS  in  the  Laplace  and  Fourier  transform  domains  that  fractional 
derivatives  are  the  easiest  to  work  with.  In  the  Laplace  transform  domain 
the  fractional  derivative  operator  with  zero  initial  conditions  has  the 
property 


LlD^Ix(t)]]  -  s'^Uxlt)]  (5) 

where  L  is  the  Laplace  transform  operator 

oo 

Llx(t)) »  |x(t)e*^^dt  (6) 

0 

Similarly,  in  the  Fourier  transform  domain. 

FlD''[i(t)l)  -  (ia))''F{x(t)l  (7) 

where  F  is  the  Fourier  transform  operator 
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00 


F|i(l)]  -  |i(t)e"i“^dt  (8) 

0 

To  gain  some  insight  into  how  the  fractional  derivative  relates  to  the 
traditional  derivative,  consider  a  sinusoidal  signal, 

w(t)  -  sinlcot)  (9) 

The  first-order  derivative  of  this  is  simply 

D*|w(t)l -wcos((«)t)-(i>sin(<«Jt  +  ^)  (10) 

So,  the  derivative  causes  a  change  in  magnitude  o  and  a  phase  shift  ( ^ 

radians  or  90*)  from  the  initial  signal.  The  fractional  derivative,  in  this  case 
the  half  order  derivative,  displays  a  similar  behavior  (asymptotically)  as  t 
approaches  infinity  and  with  zero  initial  conditions 

i  jj 

lim  D2  [w(t)]  -  sin(<i>t  +  7)  (11) 

t->»  ^ 

The  restriction  of  t  approaching  infinity  (steady-state  motion)  occurs  because 
of  transients  introduced  with  the  fraaional  derivative.  The  half  m'der 
derivative  seems,  in  a  certain  sense,  to  lie  "between”  the  signal  and  its  first- 
order  derivative  (See  Figure  1 ).  Its  phase  shift  is  45*.  midway  between  the 
unshifted  original  signal  and  the  90*  phase  shift  of  the  first-order  derivative. 
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Its  magnitude  is  scaled  by  ya .  compared  with  u  for  the  first-order 
derivative.  It  looks  as  if  the  half-order  derivative  could  almost  be 
interpolated  between  these  two  signals.  Heuristically  this  is  comforting,  since 
it  demonstrates  that  the  behavior  of  the  fractional  derivative  bears  some 
relation  to  the  more  familiar  mathematics  of  integer  order  calculus.  This  is 
not  a  completely  arbitrary  example  since  by  using  a  Fourier  series 
expansion,  many  piecewise  continuous  periodic  functions  can  be  written  as  a 
sum  of  sinusoids. 


Figure  1 .  Comparison  of  Fractional  And  Integer  Order  Differentiation 
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Overview 

In  order  to  implement  the  fractional  calculus  control  scheme,  it  was 
necessary  to  build  some  type  of  circuit  that  was  capable  of  fractional 
differentiation  or  integration.  Because  of  problems  with  low-frequency  drift 
present  in  integration  circuits,  it  was  decided  to  concentrate  on  building 
differentiators. 

Two  circuit  designs  were  selected  from  the  open  literature,  and  both 
types  were  built  and  tested  to  see  which  had  better  performance.  The 
minimal  criteria  for  success  were  chosen  to  be  a  phase  shift  of  45"'  ±5’.  and  a 
slope  of  10  ±  0.5  dB/decade  on  a  log-magnitude  plot,  within  a  frequency 
range  of  1  to  1 00  hertz. 

Since  accelerometers  were  chosen  as  the  sensors  in  this  experiment,  it 

was  necessary  to  first  integrate  the  acceleration  signal  twice  to  obtain 

position  and  velocity  measurements  of  the  response.  Integer  order 

integrators  with  good  low  frequency  response  were  available  from  a 

i  1 

previous  experiment  (9).  so  this  was  not  a  problem.  The  D2  (z(t)l  and  D2  |z(t)] 
signals  were  obtained  by  passing  the  position  and  velocity  signals, 
respectively,  through  half-order  differentiation  circuits. 
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Oldham-Zoski  Circuit 

The  first  circuit  that  was  examined  was  developed  by  Oldham  and  Zoski 
(21 )  for  use  in  electrochemical  analysis.  The  basic  component  of  the  circuit  is 
the  domino  ladder  network,  which  consists  of  interconnected  chains  of 
resistors  and  capacitors  (See  Figure  2). 


Figure  2.  Domino  Ladder  Qrcuit 


The  magnitude  of  the  potential.  Ej(t),  generated  across  the  jth  pair  of 
components,  Cj  and  Rj,  when  a  current,  i(t),  is  passed  through  the  entire 
network  is  given  by 


i(t)  - 


Ri 


(12) 
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Eq  ( 1 2)  may  be  solved  to  find 


Ej(t)  -  p 

Lj 


i(t-T)  exp 


RjCj 


dt 


(13) 


where  i  is  an  variable  of  integration.  Since  the  current,  i,  is  common  to  all 
pairs  of  components,  it  follows  that  the  potential  across  the  entire  network  is 


(14) 


At  this  point,  adjacent  resistors  and  capacitors  are  taken  to  differ  in  value  by 
respective  factors,  g"^  and  G‘j.  The  values  of  the  resistors  and  capacitors  form 
respective  geometric  progressions. 


Rj-g'jRo  (15) 

and 

Cj-G-lCo  (16) 

Here  g  and  G  are  the  geometric  factors,  both  of  which  are  taken  to  be  greater 
than  unity.  Eq  ( 14)  may  be  restated  (21:  29)  as 
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where  v  is  defined  as 

In  G 
’ln(Gg) 


(18) 


This  definition  limits  v  to  less  than  unity.  In  the  limit  as  n  and  N  both 
approach  infinity  and  g  and  G  both  approach  unity,  Eq  ( 17)  may  be  written 
as  (21:  29) 


t 

n  csc(vn)  f 

E(t)  - ; — 

r(l-v)  ln(Gg)  Co'-''  J 

0 


n  csc(va)  Rq 
In(Gg)  Co*-''  dt''-* 


i(t)-  Ew(t) 


(19) 


where  the  second  equality  is  a  consequence  of  the  Riemann-Liouville 
definition  of  the  fractional  integral  (21:  48).  The  ideal  output  voltage,  Ei^,  is 
seen  to  be  proportional  to  the  (v- 1 )  differintegral  of  the  input  current.  From 
the  principle  of  continuity,  the  input  current  equals  the  output  current,  so  if 


10 


a  resistor  is  connected  in  series  with  this  circuit,  the  voltage  across  the 
resistor.  Eult),  is  related  to  the  input  voltage  E(t)  by  (10) 


ER(t)  -  i(t)  R  - 


In(Gg)  Cq*-'' 
ji  csc(vn)  Ro'^ 


d^E(t) 

dt'' 


(20) 


Because  of  the  limit  on  v,  the  order  of  differentiation  is  limited  to  between  0 
and  1.  The  exact  order  of  the  derivative  is  determined  by  the  relationship 
between  G  and  g.  which  do  not  necessarily  approach  unity  at  the  same  rate. 

N  and  n  were  stipulated  to  be  infinite  in  the  derivation  of  Eq  ( 19).  but 
this  requires  using  an  infinite  number  of  resistors  and  capacitors,  which  is 
not  a  realistic  requirement  when  physically  constructing  a  circuit.  The 
requirement  that  g  and  G  approach  unity  for  each  resistor-capacitor  pair  also 
results  in  an  infinite  number  of  components  (2 1;  3 1 ).  So  that  a  circuit  may  be 
construaed,  n  and  N  are  chosen  to  be  much  less  than  infinity,  specifically 
n-N-2  for  the  case  of  v  -  0.5,  g  -  G  -  2.57  (21:  32).These  approximations 
result  in  upper  and  lower  frequency  bounds  on  the  circuit,  and  cause  some 
minor  distortions  within  the  bandpass  of  the  circuit.  However,  within  the 
bandpass  the  circuit  still  provides  a  reasonable  approximation  of  a  fractional 
order  derivative. 

The  actual  design  of  the  circuit  is  somewhat  complicated  by  the  addition 
of  an  extra  resistor -capacitor  pair  at  each  end  of  the  circuit.  These  are  added 
to  help  compensate  for  the  infinite  number  of  components  that  were 
truncated  when  finite  n  and  N  were  selected,  and  much  improve  the 
performance  d*  the  circuit.  For  complete  details,  consult  (21).  Detailed 
information  on  the  construction  of  the  circuits  is  contained  in  Appendix  A. 
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Qldfieid-Burd-Doe  Circuit 


The  second  circuit  eiamined  was  developed  by  Oldfield,  Burd  and  Doe 
( 19)  for  use  in  experimental  heat  transfer  work.  It  is  conceptually  very 
similar  to  the  Oldham-Zoski  design,  but  has  the  advantage  of  considering  real 
world  component  characteristics.  Specifically,  since  it  is  much  more  difficult 
and  expensive  to  obtain  high-precision  capacitors  than  resistors,  this  design 
makes  allowances  for  capacitors  that  are  only  of  minimal  precision.  The 
Oldfield-Burd-Doe  (henceforth  called  OBD)  design  was  done  specifically  to 
solve  for  the  heat  transfer  rate  in  shock  tube  experiments.  Specifically,  the 
equations 


and 


dx 


q-  -k  — 
^  ax 


(21) 

(22) 


are  solved  for  d  in  terms  of  T  (in  the  Laplace  domain)  to  yield 

1  i 

q(s)  -  (pck)2  T(s)  s2  (23) 


where  the  q(s)  and  T(s)  are  Laplace  transformed  variables.  From  the 
discussion  in  section  II.  it  is  obvious  that  this  equation  has  a  solution 
containing  a  half  order  derivative.  Since  the  designers  of  this  circuit  were 
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only  interested  in  this  one  specific  problem,  the  original  development  did  not 
contain  a  parameter  to  allow  for  the  design  of  an  arbitrary  order  fractional 
derivative.  Since  this  experiment  chose  a  priori  to  use  half  order  derivatives, 
this  did  not  present  a  problem. 

Whereas  Oldham  and  Zoski  allowed  the  progression  of  the  resistors  and 
capacitors  in  their  design  to  vary  independently,  Oldfield,  Burd  and  Doe 
designed  a  circuit  where  both  progressions  varied  at  the  same  rate,  so  g  -  G 
in  this  development.  In  exactly  the  same  manner  as  £q  ( 18)  this  limits  the 
order  of  the  derivative  to 


In  g  In  g  X 
ln(g2)  2  In  g  2 


(24) 


The  OBD  circuit  begins  with  the  same  assumption  as  before,  that  is.  an 
infinite  number  of  components  continuously  varying  in  resistance  and 
capacitance.  The  basic  circuit  is  arranged  somewhat  differently  (Figure  3). 
but  the  result  is  the  same. 


The  complex  impedance  across  such  an  infinite  network  is  given  by  ( 19:  237) 
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(25) 


z(s)  -  -yjl  A 

Once  again,  the  assumptions  necessary  for  this  result  are  relaxed  to 
allow  a  finite  number  of  components.  Rather  than  a  choosing  an  arbitrary 
geometric  progression  as  Oldham  and  Zoski  do,  Oldfield  chooses  a  progression 
for  the  capacitors  so  as  to  match  commercially  available  components.  One  of 
the  standard  series  of  capacitors  commercially  available  is  (1.0,  2.2,  4.7,  10.0, 

i 

22.0, ...).  By  choosing  the  correct  value  for  g  (103),  this  progression  can  be 
approximated  closely. 

The  Oldfield  paper  (19)  showed  experimentally  that  precise 
approximation  of  the  infinite  analog  was  not  essential,  as  long  as  the  correa 
relationship  was  maintained  between  adjacent  resistor /capacitor  pairs.  Thus, 
rather  than  selecting  a  specific  capacitance  and  trying  to  match  it  with 
various  combinations  of  relatively  inaccurate  capacitors  (21),  they 
recommend  selecting  a  single  capacitor  that  is  in  the  vicinity  of  the  desired 
value,  and  then  calculating  the  desired  resistance  to  match  it.  Since  resistors 
of  five  percent  and  even  one  percent  aoiuracy  are  available  commercially  at 
relatively  low  cost,  this  is  much  easier  and  more  inexpensive  than  trying  to 
match  both  resistors  and  capacitors  to  predetermined  values  as  the  Oldham- 
Zoski  circuit  requires. 

The  OBD  circuit  requires  some  active  components  because  the  output  of 
the  network  is  the  half  order  derivative  of  the  output  current  rather  than 
the  output  voltage  as  desired.  From  Ohm's  law  and  Eq  (23)  it  can  follows  that 
the  output  current  of  the  circuit  is  related  to  the  input  voltage  by 


M 


I(s)  = 


(26) 


■\[^E(S)  A 

A  simple  operational  amplifier  circuit  known  as  a  current-to-voltage 
converter  (26:  439)  transforms  this  into  the  desired  voltage  signal. 

The  current-to-voltage  converter  has  an  inverting  configuration  and 
causes  the  signal  to  be  180’  out  of  phase.  Also,  the  gain  of  the  circuit  is  not  as 
large  as  is  desired.  To  match  the  theoretical  fractional  differentiator,  the 
circuit  should  have  a  unity  gain  for  an  input  signal  with  an  angular 
frequency  of  one  radian  per  second.  Placing  an  inverting  amplifier  (26;  436) 
as  the  final  component  of  the  circuit  corrects  the  phase  shift  and  amplifies 
the  output  so  that  the  circuit  has  the  desired  gain  characteristics.  The  final 
configuration  of  the  OBD  circuit  is  shown  in  Figure  4.  Details  of  the  OBD 
construction  can  be  found  in  Appendix  B. 


15 
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Figure  4.  Schematic  of  OBD  Fractir  lal  Differentiator 


Characterization  of  Circuits 

In  order  to  decide  which  of  these  two  circuit  designs  was  superior,  it 
was  necessary  to  experimentally  test  the  fabricated  circuits  and  compare 
them  to  the  ideal  half -order  differentiator.  The  ideal  half -order 
differentiator  has  a  43*  phase  shift  and  10  dB/decade  slope  on  a  Bode  plot. 
To  obtain  a  transfer  function  for  each  circuit,  they  were  measured  by  a  B&K 
2032  Signal  Analyzer  with  a  band-limited  pseudorandom  noise  input. 
Typical  Bode  plots  for  the  half  order  Oldham-Zoski  and  OBD  circuits  are 
shown  in  Figures  3  and  6.  Using  similar  components  and  care  in  construction, 
the  OBD  circuit  was  found  to  be  more  accurate  than  the  Oldham-Zoski  circuit. 
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and  accurate  over  a  wider  frequency  range.  All  further  fractional 
differentiators  were  constructed  according  the  OBD  method. 
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Figure  6.  Frequency  Response  Func±ion  of  OBD  Circuit 


IV.  Characterization  of  Eiperimental  EQuioment 


Overview 

The  structure  studied  in  this  experiment  is  a  modification  of  the 
Advanced  Beam  Experiment  originally  developed  by  the  Structural 
Dynamics  branch  of  Wright  Laboratory  to  study  some  of  the  problems 
associated  with  large  space  structures  (13).  It  has  low  frequency,  lightly 
damped  modes  of  vibration  that  are  closely  spaced.  In  its  original 
configuration,  all  components  were  selected  with  an  eye  towards  eventually 
being  implemented  in  a  space  environment,  so  proof-mass  actuators  were 
used  to  control  the  structure.  Because  of  the  numerous  technical  difficulties 
involved  with  these,  they  were  replaced  in  this  experiment  by  an 
electrodynamic  shaker,  which  displayed  much  more  regular  behavior.  Since 
the  point  of  this  particular  research  was  to  examine  the  effectiveness  of 
fractional  order  feedback,  and  not  to  attempt  to  solve  all  the  problems 
inherent  in  the  control  of  large  space  structures,  this  was  not  a  grave 
liability. 

Structure 

The  structure  used  in  this  experiment  is  simply  a  six  foot  inverted 
cantilever  beam  with  a  1"  x  3/4“  rectangular  cross  section.  For  this  study  the 
circular  disc  originally  mounted  at  the  base  of  the  cantilever  beam  has  been 
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removed.  The  structure  is  instrumented  vith  two  piezoresistive 
accelerometers  and  controlled  by  a  single  electrodynamic  shaker.  The  beam 
is  arranged  as  shown  in  Figure  7.  The  physical  characteristics  of  the  beam 
are  shown  in  Table  1. 


Table  1.  Struaure  Physical  Properties 


Property  Description 
Beam  Length  (L) 

Y  Cross-Section  Width  (a) 
Z  Cross-Section  Width  (b) 
Cross-Section  Area  (A) 
Young's  Modulus  (E) 
Shear  Modulus  (G) 

Beam  Density  (p) 

Beam  Mass  (mi,) 

Y  Moment  ol  Inertia  (ly) 
Z  Moment  of  Inertia  (I,) 


Value 

Units 

70.75 

in 

1.01 

in 

0.758 

in 

0.7656 

in2 

10.8  X  10^ 

lbf/in2 

4.1 1  10^ 

lbf/in2 

2.591 1  10"* 

lbf>sec2/in^ 

1.4031  10-2 

Ibf-sec2/in 

3.6671  10-2 

in^ 

6.5081  10-2 

in^ 
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Figure  18.  Diagram  of  Experimental  Structure 
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Because  of  the  low  frequencies  of  interest  for  this  structure, 
piezoresistive  accelerometers  were  chosen  to  measure  the  response  of  the 
beam.  Piezoelectric  accelerometers  are  more  commonly  used  in  vibration 
studies,  but  are  not  as  effective  at  such  low  frequencies.  The  specific 
accelerometers  used  in  this  eiperiment  were  Endevco  Model  2262 
piezoresistive  accelerometers  with  Endevco  Model  4423  signal  conditioners. 
These  provide  good  low  frequency  response,  and  can  actually  make  static 
(DC)  measurements.  The  signal  conditioners  provide  a  zero  adjustment  so 
that  the  output  of  the  accelerometers  can  be  nulled  before  testing.  For  more 
details  see  (10)  and  (11). 


Force  measurements  were  made  with  a  B&k  Model  8200  force  gauge. 
Force  measurements  were  used  in  the  experimental  modal  analysis 
procedure,  and  taken  off  of  an  impact  hammer  used  to  excite  the  structure. 
When  using  an  actuator  to  control  a  system  it  is  common  practice  to  measure 
the  output  of  the  actuator  directly  and  feed  that  information  back  into  the 
controller.  This  is  done  because  many  actuators  display  some  sort  of 
nonlinear  behavior  at  some  point  in  their  response,  and  there  may  be  a 
substantial  difference  between  the  commanded  input  and  the  input  the 
actuator  actually  applies  to  the  system.  In  order  to  best  control  the  system,  it 
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is  necessary  to  know  precisely  what  forces  are  acting  on  the  system  to  effect 
its  reponse.  Because  we  have  only  one  force  sensor  available  for  this 
experiment,  the  actuator  was  calibrated  with  the  force  sensor,  but  run  open- 
loop  in  the  tests.  The  actuator  used  is  very  nearly  linear  in  the  frequency 
bandwidth  of  interest  so  this  did  not  prove  to  be  a  problem  for  this 
experiment. 


Actuation 

The  only  actuator  used  in  this  experiment  is  an  APS  Model  1 13-LA 
Structural  Dynamics  Shaker,  driven  by  an  APS  Model  1 1 4  dual  mode  power 
amplifier.  This  shaker  is  ideal  for  work  in  low  frequency  vibration  studies 
because  of  its  long  stroke  limit  of  six  inches  and  its  light  weight  armature, 
which  minimizes  its  influence  on  lightweight  structures.  The  shaker’s  force 
output  is  linear  in  the  frequency  range  of  this  system  (Figure  8),  so  it  was 
possible  to  run  it  open-loop,  with  no  fwce  measurement  fed  back  from  the 
attachment  point.  For  more  details,  see  ( 1)  and  (2). 

Digital  Controller 

The  controller  for  this  experiment  is  the  PC- 1000  Systolic  Array 
Processor,  which  is  a  sixteen  channel  digital  controller.  It  has  internal 
analog-to-digital  (A-D)  and  digital-to-analog  (D-A)  converters,  and  allows 
input  and  output  gains  to  be  set  for  each  channel.  It  has  a  maximum 
sampling  rate  of  2000  Hz,  which  was  the  rate  chosen  for  this  experiment. 
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The  controller  is  programmed  via  a  GPIB  connection  to  a  PC.  For  more 
details  see  (27). 


Modsling 

One  of  the  fundamental  difficulties  in  designing  any  control  system  is 
developing  a  mathematical  model  that  adequately  represents  the  physical 
system  to  be  controlled,  and  this  problem  proved  to  be  no  exception.  For  the 
purposes  of  this  experiment  only  bending  motion  in  one  axis  was  studied, 
and  it  was  assumed  to  be  uncoupled  to  bending  motion  in  the  other  axis  or  to 
torsional  motion.  original  intention  was  to  use  a  22  degree-of -freedom 
finite  element  ^del  that  had  been  previously  developed  for  the  structure 
by  Yang  (50).  This  model  looked  at  displacements  and  rotations  at  eleven 
point3  on  the  beam,  with  point  masses  representing  the  accelerometers  and 
the  shaker  armature.  The  model  was  iteratively  tuned  by  varying  a  linear 
and  torsional  spring  constant  at  its  base  until  it  produced  good  agreement  for 
the  experimentally  measured  natural  frequencies  of  the  first  two  modes. 
Details  of  the  model  are  shown  in  Appendix  F.  However,  controllers  that 
were  designed  using  this  model  had  very  poor  performance  when  actually 
tested  on  the  struaure. 

A  closer  comparison  between  the  model  and  experimental  data  showed 
that  the  model  was  only  accurate  for  the  first  two  modes,  and  of  varying 
accuracy  for  the  other  twenty  modes.  In  this  case  the  LQ  regulator  (designed 
in  spatial  variables)  weighted  all  modes  equally  when  looking  for  a  solution. 
The  resulting  controllers  had  a  very  wide  bandwidth,  which  tends  to  make  a 
system  very  susceptible  to  any  noise  present  in  the  feedback. 
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Since  only  the  first  two  bending  modes  were  to  be  controlled,  it  was 
decided  to  reduce  the  model  to  reflect  this.  This  model  would  only  look  at 
the  two  positions  where  accelerometers  were  present.  The  original  stiffness 
matrii,  ko,  taken  from  the  finite-element  model,  was  reduced  by  first 
inverting  it  to  find  the  compliance  matrii,  then  picking  out  the  four  terms 
that  affected  the  X5  and  Xi,  positions.  This  2  by  2  matrix  was  then  inverted  to 
find  the  reduced  stiffness  matrix,  k. 


k 


■  263.87  -54.82  ‘ 
.  -54.82  17.35  . 


(Ib/inch) 


(27) 


To  produce  a  better  representation  of  the  actual  behavior  of  the 
structure,  it  was  decided  to  calculate  a  reduced  mass  matrix,  m,  and  a 
reduced  damping  matrix,  c,  using  the  first  two  mode  shapes,  the  natural 
frequencies,  and  the  resonant  response  determined  experimentally.  The 
mode  shapes  are  calculated  from  the  imaginary  portion  of  the  transfer 
function  (12:  195-198).  Details  are  shown  in  Appendix  F. 

To  begin  with,  it  is  necessary  to  show  the  relationship  between  the 
displacement  formulation  and  the  modal  formulation  of  the  problem.  The 
displacement  formulation  is  based  on  a  physically  measurable  displacement 
vector,  £.  For  a  damped  N  degree-of-freedom  system,  the  motion  is  described 
by 


mi^  ci  *  Du 


(28) 
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For  this  problem,  u,  the  force  input  to  the  structure,  is  a  scalar  quantity.  D 
shows  at  which  coordinate  the  control  force  is  being  applied. 

-  1  1  0  1  (29) 

There  is  an  equivalent  formulation  written  in  modal  coordinates,  which 
are  related  to  the  physical  displacements  by 

(30) 

where  ^  is  the  vector  of  modal  displacements  and  •  is  the  modal  matrii 
composed  of  the  eigenvectors  of  the  open  loop  system.  The  equivalent 
description  is  written  as 

III  ♦l2CnWftlh+  (Wo^ln  -  •'^Du  (31 ) 

where  [Z^oWfi]  and  [wq2]  are  diagonal  matrices  containing  the  associated 
values  for  each  mode.  This  formulation  presupposes  that  the  modal  matrii, 
has  been  scaled  so  that 

(32) 

Comparing  Eq  (28)  and  Eq  (3 1 )  shows  that  k  must  be  equal  to 

♦T  k  ♦  -  (wo21  (33) 
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So  having  determined  9  and  (wq^}  experimentally,  and  with  the  t  calculated 
in  Eq  (27)  it  becomes  possible  to  scale  the  ♦  so  that  Eq  (33)  is  approximately 
true.  There  are  only  two  parameters  to  vary,  the  scaling  factors  for  each 
mode  shape,  so  it  is  in  general  not  possible  to  force  the  diagonal  terms  to 
zero  in  Eq  (33).  It  is  possible  to  scale  •  such  that  the  diagonal  terms  are  met 
exactly.  By  rearranging  Eq  (32)  it  is  now  possible  to  calculate  an  m  that 
satisfies  the  experimentally  measured  data. 


m  -  (♦T)-! 


9.28  1.83  ■ 
.  1.83  2.74  . 


X  10"^ 


(  Ibf  s2  \ 


(34) 


Similarly,  it  is  possible  to  calculate  a  damping  matrix,  c,  using  the  scaled  ^ 
and  the  diagonal  matrix  (2^aWc! 


C-  (♦T)-!  |2CnWn] 


3.46 

1.69 


1.69 
1.66  J 


X  10- 


Ibf  s  ) 


(35) 


In  order  to  use  the  tools  available  from  modern  control  theory,  it  is 
necessary  to  put  Eq  (28)  into  the  standard  state  space  form  of 


i  -  Ai  +  Bu 


(36) 


where  i  is  defined  as  the  state  vector  composed  d*  i  and  2.  and  A  and  B  have 
the  form 
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0  I 
[  -iii-*k  -m-*c 


0  0  10 
0  0  0  1 
-37300  8240  -8.22  0.606 

44800  -11800  3.56  -5.95 


(37) 


(o  0  m-iDl  -  [  0  0  124  -82.7  ) 


(38) 


The  product  -m'*k  was  calculated  from  the  definitions  of  m  and  k  in 
Eqs(32)and(33). 


-m'*k  -  -  (4>T)-1  <^l[wn2J  -  -  ♦  [Wn^J  ♦'I  (39) 

To  verify  the  accuracy  of  the  model,  the  open-loop  transfer  function 
computed  from  this  model  is  plotted  against  the  experimentally  determined 
open-loop  transfer  function  in  Figure  9.  All  transfer  functions  derived  from 

the  experimental  structure  shown  in  this  work  will  be  accelerances 

with  the  force  measurement  taking  place  at  the  midpoint  of  the  beam  where 
the  actuator  is  located  (x^),  and  the  acceleration  measurement  taking  place  at 
the  tip  accelerometer  (Xn).  The  models  show  reasonable  agreement  in  terms 
of  the  both  the  frequencies  and  dampings  for  both  modes,  but  the 
magnitudes  are  off  by  a  significant  amount. 

The  most  likely  source  of  the  magnitude  error  is  the  poor  accuracy  of 
the  measurement  technique  at  low  frequency.  The  structure  has  bending 
modes  in  the  other  lateral  direction  and  torsional  modes  as  well,  all  within 
the  bandwidth  of  inte.  .st.  The  frame  supporting  the  structure  has  modes 
occuring  in  that  frequency  range  as  well.  The  theoretical  model 


z(s)^ 

F(s)| 
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(gp)  apn-»ra3Bpi 


31 


Figure  9, 


displays  the  expected  behavior  of  a  two  degree  of  freedom  system,  with  the 
slope  approaching  40dB/decade  at  low  frequencies  and  the  measurement 
has  been  checked  thoroughly.  The  logical  conclusion  is  that  the  structure’s 
response  is  not  well  represented  by  a  second  order  mass  and  stiffness  model. 
To  better  represent  the  behavior  of  the  structure  within  the  desired 
frequency  range,  higher  order  models  should  be  examined. 
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V.  Measurement  Channels 


Overview 

A  major  source  of  difficulty  in  performing  this  experiment  was  the 
construction  and  troubleshooting  of  the  measurement  channels.  The  signal 
from  the  beam  is  acquired  as  an  acceleration,  then  it  is  variously  integrated, 
differentiated,  and  filtered  to  obtain  the  requisite  position,  velocity,  and 
fractional  derivative  feedback  signals.  A  diagram  of  the  entire  process  is 
shown  in  Figure  10. 

Integration  Circuits 

The  integration  circuits  used  in  this  experiment  were  originally  built  by 
Jacques  ( 13)  for  his  thesis  work,  based  on  a  design  provided  by  the 
Structural  Dynamics  Branch,  Wright  Laboratory  (shown  in  Figure  1 1 ).  This 
design  provides  good  performance  for  integration  above  1  Hz,  and  attenuates 
inputs  below  0.1  Hz.  where  the  pole  of  the  integrator  is  located. 

High-Pass  Filters 

To  obtain  the  position  measurements  required  for  full-state  feedback, 
the  output  of  the  accelerometers  had  to  be  passed  through  two  successive 
integration  processes.  Theoretically  there  is  no  reason  why  this  should 
present  a  problem,  nevertheless,  when  non-ideal  components  were  used 
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Figure  10.  Measurement  Channel  Diagram 


there  was  a  complication.  Real  integrators  are  normally  construaed  with  a 
transfer  function  of  the  form  ^  ~  ^  ,  as  opposed  to  that  of  the  ideal 

integrator,  which  is  7.  When  two  integrators  were  cascaded  (the  output  of 

one  fed  into  the  other),  the  transfer  function  for  the  system  is  of  the  form 


G(s) 


1 

s^  +  2bs  + 


(35) 


From  the  theory  developed  to  study  second-order  systems,  this  corresponds 
to 


35 


1 


(41) 


G(s)  -*7 

+  24WnS  +  Wq 


Comparing  these  two  expressions  reveals  that  b  is  equal  to  the  natural 
frequency  (w„)  of  the  system  and  that  the  system  is  critically  damped 
(C  -  1 ).  The  settling  time  for  such  a  system.  T,.  is  defined  as  the  time 
required  for  the  system  to  settle  within  a  certain  percentage  of  the 
amplitude  of  an  input  step  function.  For  a  2%  settling  time  this  is 
approximately  (18: 236-237) 


T. 


4  4 

^w„  '  (l)(.l) 


*  40  seconds 


(42) 


Since  the  fastest  mode  of  vibration  of  interest  in  this  structure  is  35  Hz.  this 
is  much  too  long  a  settling  time  for  the  position  measurement  taken  to  be  of 
any  possible  use.  In  fact,  the  inherent  damping  present  in  the  beam  would 
cause  it  to  settle  before  an  initial  measurement  of  its  motion  would! 

The  solution  to  this  problem  lay  in  a  closer  examination  of  the 
phenomena  involved.  The  resonance  in  the  first  integrator  is  feeding  directly 
into  the  second  integrator  and  causing  a  “ringing"  in  the  circuit.  To  avoid  this, 
a  high-pass  filter  is  placed  between  the  first  integrator  and  the  second  to 
filter  out  the  response  in  the  vicinity  of  0.1  Hz. 

Because  of  the  rapid  roll-off  (the  rate  of  attenuation  of  a  signal  with  a 
frequency  outside  of  the  pass  band)  required  of  the  filter,  it  was  not 
practical  to  build  a  passive  filter  for  this  experiment.  There  are  many 
possible  choices  of  active  filter  designs  available,  including  elliptic,  parabolic. 
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Bessel,  Papoulis,  Gaussian,  and  Butterworth  filters  (24),  all  of  which  were 
considered  at  one  time  or  another. 

Unfortunately,  the  filter  creates  another  problem.  Active  electronic 
components  do  not  function  well  as  their  inputs  approach  direct  current 
(0  Hz).  Also,  in  order  to  minimize  the  phase  distortion  of  the  input  signals  the 
corner  frequency  of  the  filter  must  be  placed  approximately  two  decades 
below  the  lowest  frequency  of  interest.  As  initially  configured,  the  structure 
had  a  first  mode  bending  frequency  of  2.1  Hz,  and  it  proved  impossible  to 
build  a  filter  that  did  not  distort  the  phase  of  the  first  mode  signals  by  an 
unacceptable  amount.  To  work  around  this,  the  structure  was  physically 
modified  by  removing  a  plate  that  had  been  attached  to  the  bottom  of  the 
beam.  This  raised  the  first  mode  frequency  to  6.3  Hz,  and  filters  were  able  to 
be  constructed  that  better  preserved  the  phase  information  at  that 
frequency. 

To  improve  the  roll-off  of  a  filter,  one  would  normally  build  a  higher 
order  filter  (add  more  poles  to  the  filter).  However,  because  of  the  extremely 
low  frequencies  involved,  this  strategy  did  not  work.  In  fact,  when  the  filter 
order  was  changed  from  second  to  third  order  (or  higher)  the  performance 
actually  deteriorated.  Since  optimal  performance  seemed  to  occur  with  the 
second  order  filters,  it  was  necessary  to  choose  a  filter  design  that  had  the 
highest  inherent  roll-off.  Chebyshev  filters  (24)  have  a  higher  roll-off  than 
the  other  filter  designs  considered,  but  cause  a  small  ripple  in  the  magnitude 
of  the  pass  band.  This  was  deemed  to  be  an  acceptable  design  trade-off. 

A  Chebyshev  filter  has  a  transfer  function  based  on  a  Chebyshev 
polynomial  appearing  in  the  denominator  of  the  transfer  function.  A  Bode 
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Figure  12.  Frequency  Response  of  High-Pass  Filter 
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Figure  13-  High  Pass  Filter  Circuit  Diagram 

plot  of  one  of  the  filters  used  is  shown  in  Figure  12.  A  circuit  diagram  of  a 
second  order  Chebyshev  filter  is  shown  Li  Figure  13.  Specific  details  of  the 
construction  are  included  in  Appendix  C. 

Ley-Pass  FUtgrs 

When  the  fractional  state  inputs  were  first  fed  back  from  the  structure 
and  the  controller  was  turned  on,  the  system  was  found  to  be  only 
marginally  stable,  with  audible  high  frequency  resonances  occuring.  When 
these  where  characterized,  it  was  found  that  the  structure  was  resonating  at 
harmonics  of  the  sampling  frequency  of  the  digital  controller.  The  problem 
was  finally  traced  to  the  high  frequency  behavior  of  the  fractional 


39 


differentiator  circuits,  which  did  not  roll-off  until  approximately  15  KHz.  (See 
Figure  14.)  The  differentiators  amplify  any  high  frequency  noise  by  the 
square  root  of  its  frequency.  Any  small  amount  of  high  frequency  noise 
entering  the  fractional  differentiators  had  its  magnitude  amplified  greatly, 
and  this  dominated  the  sample  taken  at  every  period. 

The  solution  to  this  probem  was  to  low-pass  filter  the  fractional  state 
channels  just  prior  to  entering  the  PC- 1000  so  that  there  was  at  least  20  dB 
of  roll-off  at  1000  Hz.  Once  again,  there  is  a  design  trade-off  between 
preserving  the  phase  information  at  the  frequencies  of  interest  and  still 
getting  the  necessary  roll-off.  For  this  problem,  second  order  Butterworth 
(20)  low-pass  filters  were  chosen  and  built.  These  filters  are  similar  to  the 
Chebyshev  high- pass  filters  discussed  earlier,  with  only  a  few  minor 
differences. 

The  Butterworth  filters  are  constructed  so  the  transfer  function  matches 
a  Butterworth  polynomial.  The  circuit  diagram  is  identical  to  Figure  13, 
except  that  the  resistors  and  capacitors  are  interchanged.  Finally,  the 
Butterworth  filter  has  a  flat  pass-band,  so  there  is  none  of  the  ripple  present 
associated  with  Chebyshev  filters.  A  Bode  plot  of  one  of  the  filters  used  is 
shown  in  Figure  1 5.  Details  of  construction  are  presented  in  Appendix  D. 

Final  Measurement  Channels 

To  maximize  the  accuracy  of  the  feedback  signals,  it  was  desirable  to 
amplify  the  signals  at  several  points  to  counteract  attenuation  caused  by  the 
integration  process.  The  PC- 1000  has  an  internal  analog-to-digital  converter, 
and  it  was  necessary  to  ensure  that  the  input  signals  at  that  point  be  near 


40 


30.00 


Log  Hag 

(dB) 


0.00 

100.00 


Phase 

(deg) 


00 


Frequency  Response  of  Low  Pass  Filter 


42 


5.  Frequency  Response  oT  Butterworth  Low-Pass  Filler 


the  minimum  full  scale  range.  To  avoid  saturating  the  integrators  and 
clipping  the  signals,  it  was  necessary  to  provide  amplificatioii  m  several 
stages.  To  avoid  the  sampling  probems  discussed  above  for  the  velocity  and 
position  feedback,  the  low-pass  filters  built  into  the  air  plITiers  for  those 
channels  were  set  to  1  KHz. 

Since  there  were  a  variety  of  components  cascaded  together  in  each 
channel,  it  was  decided  to  calibrate  the  entire  feedback  channel  at  once. 
Block  diagrams  for  the  four  different  signals  ar^  shown  in  Figure  10.  The 
circuits  were  calibrated  from  the  point  where  the  signal  from  the 
accelerometers  leaves  the  signal  conditioners  to  the  where  the  signal  was 
finally  input  to  the  PC- 1 000. 

The  calibration  was  performed  by  using  the  B&K  2032  and  its  random 

signal  generator  to  generate  the  plots  of  the  frequency  response  for  each 

channel  The  theoretical  magnitude  of  each  channel  can  be  calculated  at  a 

frequency  by  Gy,  -  ~  ,  where  v  is  the  order  of  integration  of  the  channel, 
w^ 

This  theoretical  value  was  calculated  for  an  average  frequency  of  the  system 
bandwidth.  This  number  was  divided  by  the  measured  value  at  the  same 
frequency,  and  multiplied  by  the  accelerometer  calibration  constants  to 
obtain  the  final  gain  for  e»ch  channel.  The  calibrated  gains  fcx*  each  channel 
are  listed  in  Table  2. 

The  final  frequency  response  plots  of  each  o(  the  types  of  measurement 

channel  are  shown  in  Figures  16-19.  The  high  pass  filter  before  the  second 

integrator  is  responsible  for  the  low  frequency  phase  disUx-tion  present  in 
1 

the  X  and  D2(i)  channels.  The  phase  distortion  at  higher  frequencies  in  the 
fractional  order  channels  caused  by  the  low  pass  filters  is  also  particularly 
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noticeable.  However,  the  magnitude  plots  all  exhibit  near  ideal  behavior  in 
the  frequency  interval  of  interest. 


Table  2.  Measurement  Channel  Gains. 


Channel  Number 

Signal 

Gain 

Units 

1 

D' (ij) 

7.67 

in  sec-*  V-* 

2 

D‘  (I2) 

7.19 

in  sec-*  V-* 

3 

I. 

.288 

inV-* 

4 

I2 

.277 

inV-* 

5 

D2  (l,) 

1.289 

in  sec  2  V-* 

6 

i 

02(12) 

2.915 

_i 

in  sec  2  V-* 

7 

3 

D  2  (l,) 

11.6 

.3 

in  SC";  2  V-* 

8 

0  2  (I2) 

31.95 

_3 

in  sec  2  V-* 
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The  original  intent  of  this  thesis  was  to  eiperimentally  prove  the 
validity  of  the  work  done  by  Yang(30).  He  showed  that  it  is  theoretically 
possible  to  transform  the  integer-order  (position  and  velocity)  controller 
calculated  by  the  Linear  Quadratic  Regulator  (LQR)  problem  into  an 
equivalent  controller  containing  the  desired  fractional  order  terms.  This 
method  was  approached  using  both  the  22  degree-of -freedom  finite  element 
model  and  the  reduced  two  degree-of -freedom  model. 


A  frequently  applied  form  of  optimal  control  theory  is  the  Linear 
Quadratic  Regulator.  A  regulator  is  defined  as  “a  feedback  controller 
designed  to  keep  a  stationary  system  within  an  acceptable  deviation  from  a 
reference  condition  using  acceptable  amounts  of  control."  (9:  167)  An  LQ 
regulator  provides  the  optimal  solution  to  a  linear  system  of  equations  (such 
as  the  state  space  model  developed  previously)  evaluated  on  a  quadratic 
performance  indei,  J.  J  is  written  as 

oo 

J  -  |!l^(t)0l(t)  +ii^(t)Rii(t))dt  (43) 

0 
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vhere  Q  is  a  positive-semidefinite  positive-definite  matrix  weighting  the 
states,  and  R  is  a  positive-definite  matrix  weighting  the  control  input. 
Because  there  is  only  a  single  actuator,  in  this  problem  u(t)  and  R  are  both 
scalar.  For  structural  dynamics,  a  convenient  choice  of  Q  is 


0- 


k  0 
0  ai 


(44) 


which  minimizes  the  kinetic  and  potential  energy  in  the  beam.  R  is  chosen 
according  to  how  much  control  authwity  is  available  in  the  actuators  and 
how  much  control  is  desired  over  the  structure.  When  R  is  made  smaller,  the 
amount  of  control  force  is  increased  which  intuitively  seems  like  it  should 
improve  performance.  However,  this  also  increases  the  bandwidth  of  the 
controller,  which  increases  the  sensitivity  of  the  system  to  any  noise  in  the 
measurements  or  errors  in  the  system  model. 

The  LQ  regulator  is  chosen  to  minimize  this  perfwmance  index  for  a 
particular  system.  It  is  calculated  from  the  solution,  S,  of  the  matrix  Riccati 
equation 


0  -  A^S  ♦  SA  -  SBR-^B^S  ♦  Q  (45) 

The  optimal  gain  for  the  LQ  regulator  is 

iI--R‘^B‘^Si--Gi  (46) 
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The  state  model  can  now  be  rewritten  as 


i- [a  -  BGIi  (47) 

A  caveat  is  necessary  at  this  point;  this  solution  is  optimal  only  in  a 
mathematical  sense  and  is  only  as  good  as  the  system  model  and 
performance  index  chosen.  It  will  only  give  what  is  asked  for  in  the 
performance  index  (with  respect  to  a  realistic  criterion  of  optimality),  and 
makes  no  consideration  of  robustness  issues  with  respect  to  noisy 
measurements. 

Fractional  Order  State  Feedback  Control  Theory 

Interest  in  fractional  order  state  feedback  control  grew  from  interest  in 
examining  the  effect  of  feedback  in  controlling  viscoelastically  damped 
structures.  Bagley  and  Torvik  (7)  showed  that  fractional  order  state 
equations  can  be  developed  where  the  state  vector  includes  fractional  order 
time  derivatives  of  a  structure’s  motion.  The  fractional  order  state  equations 
can  be  written  as 

D^y-Ay-BGy.  (48) 

where  y  is  the  fractional  state  vector. 

Yang  (30)  developed  a  method  of  calculating  meaningful  feedback  gains 
for  such  a  fractional  order  formulation.  This  section  summarizes  his  work.  He 
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developed  a  relation,  between  the  integer  order  and  fractional  order  state 
vectors, 


(49) 

where  #  is  a  state  transformation  matrix  which  can  be  determined  from  the 
eigenstructure  of  the  equations  of  motion  posed  in  expanded  fcx-m. 

Looking  at  the  N  degree-of-freedom  vibration  problem  described 
described  earlier,  with  damping  set  to  zero, 

B2+k2-Dit  (28) 

and  applying  the  composition  property  of  the  fractional  derivative  operator. 
Eq  (28)  can  be  transformed  into 

D“  *  D*'**  *  O*'**  *  ...  •  D*]  z(t)  (50) 

-  .  d2-»  •  _  ♦  D“]  Z(t) 

*  k  D**  2(t)  -  Du(t) 


This  equation  represents  the  same  system  in  (28),  but  expressed  in  terms  of 
fractional  order  derivatives.  To  ensure  that  velocity  terms  appear  in  the 
fractional  order  state  equations,  a  for  this  problem  is  chosen  to  be  of  the 
form  ~  where  n  is  an  integer.  This  allows  initial  value  problems  to  be  solved 

where  an  initial  velocity  is  given. 
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In  (M-der  to  solve  this  system  of  equations,  it  is  possible  to  expand  these 
N  equations  into  a  new  system  of  2  N  n  equations  of  the  form 

D«  M  i(t)  ^  l2(t)  -  F(t)  (51) 


where 


0  0  ...  0  B 

0  0  ...  B  0 


If 


0  B  ...  0  0 

B  0  ...  0  0  - 


‘  0  0  ...  0  -B  0' 

0  0  ..  -B  0  0 

0  -B  ...  0  0  0 

-B  0  0  0  0 

-  0  0  ...  0  0  k- 


"  D2-«z(t) 
02-202(1) 


i(t)- 


< 


> 

D2“i(t) 


I  l(t)  J 


(52) 


(33) 


(54) 
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0 

0 


F(t)- 


0 

0 


LDudM 


(55) 


This  augmented  set  of  2  N  n  equations  is  called  the  expanded  equations 
motion  the  system.  M  and  C  are  known  as  the  pseudo-mass  and  the 
pseudo-stiffness  matrices.  Notice  that  the  lowest  partition  of  Eq  (51 )  above 
identical  to  the  original  equation,  Eq  (28). 

Setting  the  right  hand  side  (tf  Eq  (3 1 )  equal  to  zero,  it  is  possible  to  solve 
for  the  resulting  homogeneous  system's  expanded  eigenvalues  and 
eigenvectors.  The  eigenvalues  will  appear  in  complex  conjugate  pairs.  Xi  and 

Xj.  The  eigenvectors  (^i  and  associated  with  the  conjugate  eigenvalue 
pairs)  will  be  used  to  construct  a  state  transformation  matrix  to  the  integer 
order  state  equation. 

From  the  expansion  theorem  of  theoretical  modal  analysis  it  is  known 
that  the  system  response  can  be  written  in  terms  of  the  normal  modes  of  the 
system  multiplied  by  a  set  of  time-dependent  generalized  coordinates, 

fl(t)  -  ♦ii(t)  (56) 


where 


SL(t)  -  column  vector  of  the  generalized  system  response 
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•  -  tbe  modal  matrii  consisting  of  the  modal  vectors  arranged 
in  a  matrii 

h(t)  -  a  column  vector  of  the  modal  coordinates  of  the  system. 


Applying  this  to  the  expanded  equations  of  motion,  it  is  possible  to  express 
the  system  response  as 


z(t)-*n(t) 


where 


*■[ 


I  i  ^2  I  ^2  I  -  i  I 


] 


Il(t)  - 


tll(t) 

n2(t) 

n2(t) 


nN(t) 
w  tiN(t)  ^ 


(57) 


(58) 


(59) 


Using  Eq  (57)  it  is  possible  to  extract  the  necessary  information  to 
construct  the  integer  state  vector  and  the  fractional  state  vector  for  the 
system.  The  stable  roots  of  the  expanded  equations  oS  motion  have  positive 
real  parts.  To  begin  with  selea  the  column  vectors  (<I»i)  (rf  Eq  (58)  that 
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correspond  to  the  eigenvectors  in  Eq  (39)  that  have  positive  real 
components.  Define  this  matrix  as  4^. 

The  integer  state  vector  can  be  written  as 

l(t)-^All(t)  (60) 


where  i(t)  is  defined  as 


l(t)- 


z(t) 

.  D‘Iz(t)l  . 


(61) 


and 


•a  ~  matrix  consisting  row  vectors  in  that  are 
associated  with  the  states  defined  in  Eq  (6 1 ) 

Similarly,  the  fractional  state  vector,  y.  can  be  written  as 

y(t)-*Bll(t)  (62) 


where 
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y2N-i(t) 

y2N(t)  J 


^  -  matrix  consisting  of  row  vectws  in  that  are 
associated  with  the  states  defined  in  Eq  (63) 

The  fractional  state  vector.  x(t).  is  a  nearly  arbitrary  subset  of  i  (t)  (the 
state  vector  of  the  expanded  equations  of  motion).  It  is  of  the  same 
dimension  as  the  integer  sute  veacMr,  and  may  contain  any  combination  of 
fractional  and  integer  order  states  from  any  combination  of  nodes. 

From  Eq  (62)  q  may  be  defined  as 

ll(t)  -  y(t)  (64) 

Which  can  be  substituted  into  Eq  (60) 

i(t)  -  [♦a  y(t)  -  ♦^(f)  (65) 

Which  is  the  desired  transformation  from  integer  state  vector  to  fractional 
state  veaw. 

Applying  this  result  to  the  integer  state  optimal  control  law  in  Eq  (47) 


♦  i(t)  -  [a  -  BG]  ♦x(t) 


(66) 


Which  can  be  rewritten  as 

i(t)  -  [a*  -  B*G*|  y(t)  (67) 


where 


A*  .^1  A« 

(68) 

B*  B 

(69) 

G*  -G^ 

(70) 

This  is  the  desired  result  of  a  set  of  fractional  order  state  equations  with 
optimal  LQR  feedback  control. 


Results 

There  were  several  problems  observed  when  the  LQR  method  was 
implemented.  The  finite  element  model  used  was  not  suitable  for  use  with 
LQR,  and  the  classical  integer  order  controllers  generated  with  it  were  not 
very  effective.  As  mentioned  before,  the  finite  element  model  was  only 
accurate  for  the  first  two  modes,  and  had  fifty  percent  error  in  frequency  for 
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mode  three,  with  higher  order  modes  being  progressively  more  poorly 
represented.  Since  LQR  weights  all  modes  equally  when  formulated  in 
position  variables,  the  poor  information  dominated  the  correct  information 
for  the  first  two  modes.  The  LQR  method  did  indeed  give  an  optimum 
solution,  but  it  was  for  a  structure  that  did  not  exist.  Because  of  all  the 
control  force  required  to  control  all  twenty-two  modes  of  the  model,  the 
actuator  reached  its  force  limits  before  modes  one  and  two  had  reached  even 
five  percent  critical  damping. 

The  reduced  order,  two  degree-of -freedom  model  gave  much  better 
performance.  The  integer  order  LQR  gains  were  meaningful,  and  produced 
the  desired  damping  levels  in  the  structure,  but  when  the  fractional  order 
state  vector  was  calculated,  the  gains  on  the  fractional  derivative  terms  were 
nearly  ignored.  For  example,  solving  the  LQR  problem  using  PC-Matlab  (22) 
for  R  -  0.1  yielded  integer  order  gains  of 

G-  17.90  -0.960  -0.471  0.0969]  (71) 

for  a  state  vector  composed  of 

(72) 

For  a  state  vector  composed  of 


z(t) 

D'(z(t)l 
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r 


i^(t)  »  1 


D2  IZi(t)l 
D>  lz,(t)l 

D2  lZ|(t)l 


L  z,(t)  J 


(73) 


the  state  transformation  matrix,  was  calculated  to  be 

0  0  0  1  “ 

-0.00360  -0.05351  0.331  3-50 
^  “  0  1  0  0 

-  0.331  -10.2  115  -451  - 


(74) 


Using  this  transformation  a  new  set  of  gains  is  calculated,  G*.  is  found  to  be 
G*  -  [  0.0286  -0.470  10.9  -39.2  ]  (75) 

Although  this  result  is  taken  from  measurements  at  only  the  midpoint,  it 

should  theoretically  provide  the  same  control  authority  as  the  controller 

using  position  and  velocity  measurements  from  both  points. 

3 

Another  possibility  is  feeding  back  the  velocity  and  D2  terms  from  both 
sensors.  This  would  eliminate  the  need  for  the  second  integration  process 
and  the  consequent  phase  distortion  caused  by  the  high-pass  filter.  The 
transformation  matrix  for  this  term  will  be  denoted  by  and  was  found  to 
be  equal  to 
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" -0.00166  -0.000973 

0.00935 

0.00380  “ 

-0.00530  -0.00467 

0.0207 

0.0211 

0 

0 

1 

0 

0 

0 

0 

1 

(76) 


and  the  associated  gain.  G*.  is 


G*  =  [ -0.00802  -0.00321  0.525  0.107  |  (77) 


where  the  associated  state  vector  is 


i:(t)- 


< 


D2  Iz(t)l 
D‘  lz(t)l  . 


(78) 


The  numerical  results  are  shown  below  in  Table  3,  while  the  closed  loop 
frequency  response  is  shown  in  Figure  20. 

Table  3-  Results  of  LQR  Method  (Percentage  of  Critical  Damping) 


Mode 

Predicted 

G 

G* 

G* 

1 

21.1% 

22.3% 

18.5% 

18.6% 

2 

10.9% 

8.5% 

1.8% 

7.5% 
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Closed  Loop  Frequency  Response  -  State  Transformation  Method 


(gp)  apn-^raSew 
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Figure  20.  Closed  Loop  Frequency  Response  of  LQR  Method 


The  results  shown  are  promising  but  certainly  not  conclusive.  As  in  the 
open  loop  case,  there  is  a  large  discrepancy  between  the  predicted  and 
measured  closed  loop  response.  The  shape  of  the  predicted  response 
corresponds  well  to  the  eiperimentally  measured  ones,  but  there  is  a  large 
offset  in  magnitude  that  needs  to  be  resolved  at  some  point. 

Comparing  the  integer  order  case,  G,  to  the  cases  with  fractional  order 
feedback  (G*  and  G*)  still  yields  some  useful  information,  however.  The  case 
where  all  measurements  have  been  taken  from  a  single  point  (collocated 
with  the  actuator),  G^,  is  only  marginally  stable,  with  audible  high  frequency 
resonances  occuring.  The  performance  is  degraded  for  the  first  mode,  but 
still  provides  a  reasonable  amount  of  damping.  Performance  at  the  second 
mode  is  abysmal,  and  even  destabilizes  the  system  somewhat.  (The  friction 
in  the  armature  of  the  actuator  provides  22%  damping  for  mode  2  when  the 
aauator  is  attached  but  unpowered,  so  the  performance  would  actually  be 
improved  at  this  frequency  by  turning  the  system  off.) 

The  case  where  signals  are  taken  from  both  sensors,  G*,  is  much  more 
promising.  The  closed  loop  response  is  very  similar  to  that  of  the  integer 
order  case,  although  there  is  still  some  degradation  in  the  performance  of  the 
system.  Because  of  the  known  errors  in  the  system  (modelling  errors  and 
phase  shifts  in  the  measurements),  it  is  not  possible  to  say  conclusively  that 
fractional  order  feedback  is  equivalent  to  integer  order  feedback.  However, 
it  possible  to  say  that  the  method  of  transforming  integer  order  gains  into  an 
equivalent  system  containing  miied  fractional  and  integer  order  gains  does 
seem  to  be  successful.  The  fact  that  control  remains  stable  in  the  face  of  such 
large  phase  shifts  in  the  measurements  even  indicates  reasonably  large 
phase  margins,  and  some  degree  of  robustness. 
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Overview 

An  alternative  approach  to  introducing  fractional  order  feedback  into 
structural  control  was  discovered  by  Leonard  ( 15)-  He  found  an  interesting 
result  that  seemed  to  allow  a  large  amount  of  flexibility  when  selecting  pole 
placement  gains  if  fractional  order  feedback  signals  were  available.  The 
theory  presented  here  summarizes  much  of  Leonard’s  work.  Briefly,  it 
calculates  a  set  of  equality  and  inequality  constraints  between  the  fractional 
and  integer  order  terms  that  is  necessary  for  equivalent  response  to  be 
obtained.  Once  again,  the  primary  goal  of  this  research  was  proof  of  concept, 
showing  that  fractional  derivative  feedback  would  work  in  a  real  system. 

Theory 

Assuming  an  nth  order  linear  dynamic  system  of  the  form 

i(t)  -  Ai(t)  +  Bii(t)  (79) 

where  i  contains  integer  order  feedback  terms  (position  and  velocity)  and 
there  exists  state  feedback 

il(t)  -  -Gi(t)  (80) 
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where 


G  -  |gi  g2  go| 

So  that  the  system  can  be  written 


(81) 


x(t)  -  |a  -  BG|l(t)  (82) 

Let  a(s)  be  an  equation  with  the  desired  pole  structure  of  the  system 

a(s)  -  -  Pi)  -  0  (83) 

1 

In  order  for  the  gain  matrix  to  move  the  poles  of  the  system  to  the  desired 
locations,  the  characteristic  equation  of  the  closed-loop  system  must  equal 

detjsl  -  A  +  BGj-a(s)  =  0  (8-1) 

In  other  words,  the  coefficients  of  the  characteristic  equation  must  equal  the 
coefficients  of  a(s). 

Now  suppose  that  the  equations  of  the  system  were  posed  in  fractional 
state  form. 


i 

D2  (iF(t)l  -  AiF(t)  +  Bu(t) 

where  if(i)  now  includes  the  half  order  fractional  states 


(85) 
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I(t) 

D2  lx(t)l 
D«  (i(t)J 

D2  Il(t)l 

Examining  the  pole  placement  problem  in  this  case  requires  a  little  more 
work.  Consider  the  transformation 

P  -  ^/s  (87) 

The  eigenvalue  problem  for  the  fractional  state  case  is  posed  as 

det  [  pi  -  A  ^  BG]  -  0  (88) 

The  fractional  state  system  has  twic^  as  many  poles  as  the  integer  state 
system  in  the  Vs  or  p  plane.  This  function  is  not  unique,  and  a  discontinuity 
must  exist  somewhere  in  the  space.  Poles  in  the  p  plane  with  phases  of 
greater  than  90*  or  less  than  -90*  cross  a  branch  cut,  and  will  map  to 
another  Riemann  sheet  when  they  are  brought  back  into  the  s-piane  (See 
Figure  2 1 ).  The  stability  criterion  in  the  p  plane  is  that  the  phase  of  a  pole 
must  lie  in  the  ranges  of 

45*<0<9O*  (89) 
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Figure  21.  Stability  Diagram  of  p  Plane 


or 


-45*>(^>-90*  (90) 

So  now  each  real  pole,  Pi.  may  be  factored  into  two  poles  in  the  p  plane. 
There  are  n  pairs  of  poles  and  2n  gains  so  it  appears  that  the  expanded 
problem  has  left  you  in  a  similar  situation  to  where  you  began,  except  that 
the  order  of  the  system  has  doubled.  If  this  were  true,  this  would  be  a 
remarkably  useless  result;  however,  it  is  only  necessary  to  control  half  of 
these  poles.  To  retain  similar  control  authority  to  the  integer  state  problem, 
it  is  only  necessary  to  place  n  poles  in  the  right  half  plane.  The  other  n  poles 
merely  need  to  be  kept  in  the  left  half  p  plane  off  of  the  principal  Riemann 
sheet  to  ensure  stability. 
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To  take  advantage  of  this,  first  compute  the  square  roots  of  the  desired 
pole  locations: 


(91) 


Form 


n 

fl(P)  -  n(P  ■  Pi) 
i-1 


(92) 


Form  the  characteristic  equation  of  the  system  in  the  p  plane 


Z(p)  -  det  |pl  -  A  -  BG]  -  0  (93) 

Since  Z  is  of  order  2n  and  n  is  only  of  order  n.  synthetic  division  is  possible. 
Performing  synthetic  division  on  Z(p)  with  n(p)  yields  'P(p)  and  R(p).  the 
remainder. 


z(p) 

fi(p) 


-’P(p) 


R(P) 

n(p) 


(94) 


'i'(p)  is  a  polynomial  of  order  n.  For  n(p)  divide  evenly  in  to  Z(p).  and  be 
equal  to  its  roots,  the  coefficients  of  R(p)  must  be  equal  to  zero.  This  creates 
n  equality  constraints  on  the  system. 

'P(p)  contains  the  n  poles  that  are  not  mapped  into  the  s  plane,  so  they 
must  lie  in  the  left  hand  P  plane.  To  ensure  their  stability,  it  is  possible  to 
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use  Routh-Hurwitz  criterion  to  create  inequality  constraints  that  provide 
bounds  on  stability. 

The  Routh-Hurwitz  criterion  (18:  363-365)  states  that  for  an  equation  of 
the  form 


aos®  +  ais®->  +  a2S®-2  +  a^-is  +  an  -  0  (95) 

to  be  stable,  all  coefficients  aj  must  be  of  the  same  sign  and  not  equal  to  zero, 
and,  the  coefficients  of  the  first  column  of  the  Routh  array  must  also  have 
the  same  sign. 

These  inequality  constraints  define  a  region  in  the  controller  space,  and 
any  solution  lying  in  this  region  should  provide  a  stable  solution.  This  means 
that  the  problem  has  an  infinite  number  of  non-unique  solutions,  all  of 
which  satisfy  stability  requirements.  From  a  control  system  design 
viewpoint,  this  is  a  very  interesting  result. 

At  this  point,  there  are  no  guidelines  as  to  how  best  use  this  design 
fleiibility.  However,  since  any  set  of  gains  that  satisfy  these  constraints 
should  work,  and  the  object  of  this  thesis  was  merely  to  prove  that 
fractional  derivative  feedback  is  possible,  this  method  was  acceptable,  if 
somewhat  awkward. 

Results 

Because  of  the  synthetic  division  in  this  method,  there  is  not  yet  a 
constructive  solution  to  this  problem.  Rather  than  solving  a  series  of  high 
order  determinants  and  a  lengthy  synthetic  division  by  hand,  the  symbolic 
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capabilities  of  Mathematical 29)  were  brought  to  bear  on  the  problem.  The 
script  used  to  calculate  the  equality  and  inequality  constraints  and  the 
constraints  themselves  are  shown  in  Appendix  G. 

Arbitrarily  choosing  damping  ratios  of  twenty  percent  for  mode  1  and 
ten  percent  on  mode  2,  the  desired  pole  locations  were  calculated  using 

Pi-(-CiWi  ±  Wj)  (96) 

Using  the  system  model  developed  in  Eqs  (37)  and  (3S)  the  following  integer 
order  gains  were  calculated  from  Eq  (85) 

Gp-[-3.16  2.34  0.566  0.115  ]  (97) 

Using  this  solution  as  a  starting  point,  the  integer  order  gains  were  reduced 

in  magnitude  while  the  four  dependent  fractional  gains  were  calculated 

using  the  equality  constraints.  The  resulting  set  of  gains  was  checked  versus 

the  inequality  constraints  to  ensure  that  the  closed  loop  system  would  still 

be  stable.  The  first  controller  developed  contained  all  eight  feedback  states 

3  i 

(two  each  of  position,  velocity,  D2,  and  D2  terms),  arranged  as  in  Eq  (86). 

G1  -  [  1.5  -2  1.08  0.886  0.4  0.02  0.0138  0.00743  j  (98) 

1 

The  second  controller  was  calculated  by  zeroing  the  position  and  D2  terms, 

3 

with  only  velocity  and  D2  terms  being  fed  back.  The  equality  constraints 
provided  the  necessary  information  to  calculate  the  other  four  gains.  It  was 
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thought  that  this  might  prove  to  be  a  useful  special  case  because  of  the 
difficulty  with  performing  the  second  integration  discussed  in  Chapter  3. 

G2-[  0.591  0.176  -0.00730  -0.00776  1  (99) 

These  three  controllers  were  tested  and  the  resulting  closed  loop  transfer 
functions  plotted  versus  the  theoretical  closed  loop  transfer  function  for  the 
system  in  Fiqure  22.  The  results  are  listed  in  Table  4. 

Table  4.  Results  of  Pole  Placement  Method 


Predicted 

Cp 

G1 

G2 

Mode  1 

20% 

21.6% 

19.3% 

20.1% 

Mode  2 

10% 

10.1% 

10.6% 

8.6% 

The  results  agree  fairly  well  with  the  predicted  values.  The  shape  of  the 
predicted  response  again  corresponds  well  to  the  experimentally  measured 
ones,  but  there  is  the  same  large  offset  in  magnitude  as  was  present  in  the 
previous  case. 

However,  the  most  important  comparison  to  be  made  here  is  that 
between  the  integer  order  controller,  Gp,  and  the  controllers  using  fractional 
states,  G1  and  G2.  The  point  of  this  work  was  to  show  that  control  using 
fractional  order  feedback  was  equivalent  to  schemes  using  the  more 
thoroughly  studied  integer  order  feedback.  Since  these  controllers  were 
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Frequency  Response  -  Pole  Placement 
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Figure  22.  Closed  Loop  Frequency  Response  of  Pole  Placement  Method 


calculated  from  the  same  model,  whatever  its  inaccuracies,  the  results 
should  be  very  similar.  These  numbers  show  that  this  is  indeed  true. 

Comparing  the  integer  order  feedback,  Gp,  to  the  G1  and  G2 
controllers.  G1  displays  an  error  of  1 1  %  and  5%  for  modes  one  and  two 
respectively,  while  G2  shows  an  error  of  6.9%  and  15%  respectively.  These 
results  are  encouraging,  particularly  in  light  of  the  known  errors  in  the 
circuitry.  Because  of  the  large  phase  shifts  present  in  the  fractional 
derivative  feedback  due  to  the  final  low  pass  filtering  process,  it  was 
anticipated  that  some  loss  of  control  authority  would  occur.  In  fact,  it  is 
encouraging  that  these  controllers  are  robust  enough  to  still  produce  stable 
control  in  the  presence  of  such  a  large  phase  lag. 

A  final  result  is  the  force  history  in  the  time  domain  of  the  system  with 
the  various  controllers  in  the  loop.  It  was  a  prior  concern  that  there  might  be 
some  undesirable  transients  in  the  force  response  because  of  the  effects  of 
the  poles  off  of  the  principal  Riemann  sheet.  Those  poles  that  are  close  to  the 
branch  cut  can  “bleed  over"  and  influence  the  system  behavior.  However, 
those  worries  seem  to  be  unfounded,  and  the  force  histories  of  the  different 
controllers  are  qualitatively  similar.  They  are  shown  in  Figure  23- 


73 


Figure  23. 
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Because  of  the  equality  constraints  of  the  pole  placement  method,  any 

case  with  just  four  non-zero  gains  should  be  unique  (if  a  solution  exists), 

since  the  problem  reduces  to  four  linear  equations  in  four  unknowns.  If  that 

is  so,  the  question  is,  is  there  any  correspondence  between  results  calculated 

by  the  LQR  method  and  the  pole  placement  method?  In  this  chapter,  gains 

were  calculated  for  both  the  integer  order  case  (position  and  velocity  only) 

3 

and  the  case  where  velocity  and  the  Dz  signals  were  fed  back  from  both 
sensors.  In  the  last  chapter,  a  relationship  between  those  two  gains  was 
calculated  in  Eq  (76).  so  it  is  possible  to  compare  the  two  results. 

Using  the  resuits  of  Eq  (97)  and  Eq  (76), 

Gst  -  Gp  ♦**[  0.585  0.152  -0.00715  -0.00784]  (100) 

which  compares  well  with  the  result  from  Eq  (99) 

G2-|o.591  0.176  -0.00730  -0.00776  ]  (101) 

Since  the  pole  placement  method  keeps  track  of  the  poles  that  lie  (tff  of  the 
principal  Riemann  sheet,  it  is  probably  the  more  accurate  of  the  two 
techniques.  Additionally,  it  is  the  more  general  of  the  two  techniques,  since  it 
allows  all  2  n  gains  to  be  assigned  while  the  state  transformation  technique 
of  the  last  chapter  only  allows  the  use  of  a  subset  containing  n  of  the  gains. 
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VIII.  Conclusions  and  Recommendations 


The  theoretical  predictions  made  by  Yang  (30)  and  Leonard  (15)  both 
performed  as  expected  within  the  bounds  of  the  experiment’s  error. 
Fractional  derivative  feedback  may  in  fart  prove  workable  in  a  real  world 
setting  in  the  future.  Based  on  the  results  of  this  work,  it  at  least  merits 
further  consideration  and  research. 

There  are  several  topics  that  need  to  be  addressed  before  any  more 
positive  recommendations  can  be  made.  The  first  is,  what  can  fractional 
derivative  feedback  accomplish  that  integer  order  feedback  cannot?  Feeding 
back  twice  as  many  inputs  to  get  the  same  control  authority  is  not  an  idea 
that  is  going  to  be  of  much  use  to  designers.  On  the  other  hand,  if  those  extra 
inputs  can  be  used  to  improve  system  performance,  particularly  to  increase 
system  robustness,  it  becomes  more  attractive  and  potentially  useful. 

Second,  developing  a  mass-stiffness  model  for  a  real  structure  is 
somewhat  unwieldy,  a  modal  model  is  in  general  easier  to  generate.  To  do 
this,  a  modal  estimator  for  fractional  derivatives  needs  to  be  developed. 
Eventually  a  stochastic  estimator  (a  (Caiman  filter)  would  be  desirable  to 
have  available.  Accelerometers  tend  to  be  accurate  but  not  precise,  and  if  the 
measurement  noise  statistics  are  available,  a  (Caiman  filter  can  greatly 
improve  system  performance.  Also,  since  modeling  large  space  structures  in 
a  gravity  field  is  so  difficult,  building  adaptive  abilities  into  a  control  system 
is  frequently  desirable,  and  many  adaptive  control  systems  begin  with  the 
Kalman  filter  as  a  building  block. 
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The  neii  step  for  this  particular  laboratory  setup  should  be  the 
construction  or  purchase  of  better  filters  for  the  fractional  derivative 
channels.  The  phase  shifts  on  those  channels  are  unacceptably  large,  and 
may  well  be  degrading  the  performance  of  the  controllers  using  fractional 
derivative  feedback.  Theoretically  four  fractional  derivative  feedbacks 
should  provide  the  same  damping  as  four  position  and  velocity  feedbacks  do. 
Currently  they  do  not,  there  is  some  degradation  of  the  performance.  It 
should  be  established  whether  or  not  this  is  the  nature  of  the  fractional 
derivative  feedback,  or  is  merely  the  result  of  poor  filters. 

Of  equal  or  greater  importance  is  the  resolution  of  the  errors  in 
modeling  this  system.  The  unmodelled  dynamics  are  significant  in  the 
bandwidth  of  the  system,  and  a  higher  order  model  should  be  examined. 
Once  the  a  better  model  is  in  hand,  this  work  needs  to  be  repeated, 
particularly  the  case  where  all  feedback  signals  are  generated  from  a  signal 
colocated  sensor  at  the  actuator.  This  is  potentially  the  most  useful  case  and 
it  needs  to  be  conclusively  established  whether  the  fractional  order  feedback 
IS  effective. 

The  pole  placement  technique  used  in  Chapter  7  is  not  a  design  method 
per  se.  and  it  becomes  very  complex  computationally  as  system  order 
increases.  Performing  the  synthetic  division  while  carrying  2  n  gain 
variables  requires  a  program  capable  of  symbolic  mathematical  operations 
(such  as  Mathematica  or  MACSYMA)  for  n  greater  than  four.  If  it  is  to  be 
used  on  more  complex  systems,  a  belter  method  needs  to  be  found  of 
calculating  the  controller  gains.  Algorithms  exist  that  constructively  find  the 
solution  to  the  integer  order  pole  placement  problem,  and  they  should  be 
examined  to  see  if  an  extension  to  the  fractional  order  case  is  possible. 
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Appendii  A:  Oldham-Zoski  Circuit  Design 


The  construction  of  the  Oldham-Zoski  circuit  closely  follows  that  of  the 
author’s  in  their  paper!  2 1 ).  The  design  values  are  those  given  for  the  v  -  ^ 

design.  For  more  discussion  of  this  design,  see  also  (M).  Table  A.l  lists  the 
design  and  experimental  values  for  the  test  circuit  built. 


Table  A.l  Design  and  Experimental  Values  of  Oldham-Zoski  Components 


Resistors 

(in  Mn) 

Resistor 

Design 

Measured 

Ro 

25A 

25.4 

Ri 

9.88 

9.86 

R2 

3.84 

3.83 

R3 

1.501 

1.494 

R^ 

0.5816 

0.5821 

R5 

0.2262 

0.2265 

Re 

0,08798 

0.08792 

R? 

0.00172 

0.00172 

R« 

0.00364 

0.00364 

Capacitors 

(injiF) 

Capacitor 

Design 

Measured 

Co 

10.0 

10.03 

Cl 

3.891 

3.89 

C2 

1.514 

1.511 

C3 

0.589 

0.587 

Cl 

0.2292 

0.2297 

C5 

0.0892 

0.0876 

Ce 

0.0347 

0.0351 

c? 

0.0270 

0.02707 
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The  OBD  circuit  was  designed  with  a  series  of  six  resistors  and 
capacitors  for  this  experiment.  The  capacitances  are  measured  and  then 
suitable  resistances  are  calculated,  so  their  are  no  design  values  for  the 
capacitors  per  se.  Again,  this  is  a  straightforward  application  of  the  method 
put  forth  in  ( 19).  Design  and  experimental  values  are  shown  in  Table  B.l. 

Table  B.l.  Design  and  Experimental  Values  of  OBD  Circuit  Components 


Resistors 

(in  kn) 

Capacitors(|i.F) 

Design 

Measured 

Measured 

Ri 

3.695 

3.711 

Cl 

0.0997 

R2 

14.00 

14.00 

Cz 

0.2262 

R3 

30.64 

30.58 

C3 

0.475 

64.79 

64.79 

C4 

1.012 

R5 

142.4 

142.1 

C5 

2.296 

Re 

305.3 

306.1 

Ce 

4.655 
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Aopendii  C:  Chebvshev  High-Pass  Filter  Design 

The  high-pass  filters  that  were  eventually  built  for  this  experiment 
were  the  product  of  as  much  trial  and  error  as  actual  design.  Active 
electronic  components  (such  as  op-amps)  do  not  function  well  at  extremely 
low  frequencies,  and  their  performance  becomes  erratic.  The  filter  that 
provided  the  best  performance  when  actually  built,  was  a  second  order 
Chebyshev  filter  as  shown  in  Figure  9.  The  operational  amplifier  used  was  an 
LM-747j,  which  is  a  general  purpose  dual  op-amp.  The  component  values 
that  were  found  to  work  best  were: 

C  -  0.40  |iF 
R1  -447.6  kn 
R2-  1.172  Mn 
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In  order  to  eliminate  the  resonances  caused  by  high  frequency  noise 
entering  the  sampling  process  of  the  A-D  board  on  the  PC- 1000  in  the 
fractional  derivative  feedback,  it  was  necessary  to  construct  four 
Butterworth  low-pass  filters  for  the  experiment.  The  circuit  diagram  is 
shown  in  Figure  D.l.  Once  again,  LM-747j  operational  amplifiers  were  used, 
and  the  component  values  were: 

R-  11.8  kn 
Cl  -  190  nF 
C2  -  95  nF 


The  finite-element  model  for  this  structure  is  drawn  directly  from  Yang 
(30;77-'^8)  and  models  the  beam  as  an  undamped  22  degree-of -freedom 
system.  It  models  the  beam’s  behavior  by  looking  at  lateral  motion  and 
rotation  at  eleven  points.  It  also  assumes  that  both  lateral  bending  in  the 
other  direction  and  torsional  motion  are  completely  uncoupled  from  bending 
in  the  z  direction.  The  mass  and  stiffness  matrices  are  constructed  from 
elemental  matrices.  Kei  and  Mei.  which  are  given  by 


“  12 

6Le 

-12Le 

6Le  " 

El 

6Le 

4Le2 

-6Le 

2Le2 

-l2Le 

-6Le 

12 

-6Le 

-  6Le 

2Le2 

-6Le 

4Le2- 

and 


Mel  “ 


"  156 

22Le 

54 

-13Le“ 

pALe 

22Le 

4Le2 

13Le 

-3Le2 

420 

54 

I3Le 

156 

-22Le 

-  -13Le 

-3Le2 

-22Le 

4Le2  - 

(F.2) 


where  E,  I.  p,  and  A  are  given  in  Table  1 .  and  Le  is  the  length  of  the  beam 
element  and  is  equal  to  7.075  in.  These  elements  are  summed  along  the  main 
diagonal  to  form  the  22  DOF  mass  and  stiffness  matrices.  The  mass  matrix  is 
further  modified  to  include  the  mass  of  the  accelerometers  at  positions  15 
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and  Z|  I  (0.0000343  lb  sec^  in*l )  and  the  mass  of  the  shaker  armature  at 
position  X5  (0.00185  lb  sec2  in‘l ).  The  model  is  tuned  to  match  the 
experimentally  determined  natural  frequencies  of  the  structure  by 
iteratively  changing  the  linear  and  torsional  spring  elements  connecting  the 
structure  to  its  support. 
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Apoendii  F:  Eioerimental  Modal  Analysis 


The  modal  data  of  the  structure  was  obtained  by  performing  an 
experimental  modal  analysis  using  a  calibrated  impact  hammer.  The 
procedure  has  three  major  steps:  Calibration.  Data  Acquisition,  and  Data 
Analysis.  The  resulting  data  was  used  to  construct  the  model  of  the  system, 
which  in  turn  was  used  to  calculate  the  controller  gains. 

Calibration 

Two  types  of  sensors  are  used  to  acquire  data  for  this  test, 
accelerometers  and  a  force  gauge.  The  accelerometers  used  are  of  the 
piezoresistive  type,  and  can  be  calibrated  versus  local  gravity  by  nulling  the 
sensor  when  its  placed  in  a  vertical  orientation,  then  turning  it  180*  and 
measuring  the  resulting  voltage.  The  voltage  is  proportional  to  a  2g 
acceleration  and  a  calibration  constant  may  thus  be  calculated.  Results  are 
shown  in  Table  G.l.  These  numbers  are  compared  to  the  manufacturer's 
specifications  as  a  reality  check. 

Table  G.  1 .  Accelerometer  Calibration  Data 


Accel. 

Accel. 

Signal 

2g  Meas. 

Exp.  Cal. 

Spec.  Cal. 

Position 

Serial  * 

Conditioner 

(V) 

Constant 

Constant 

Gain 

(V/g) 

(V/g) 

15 

AH72 

50 

1.764 

0.882 

0.8865 

III 

BA26 

50 

1.996 

0.998 

0.9915 
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These  numbers  compare  veil,  with  both  being  less  than  one  percent  off  of 


manufacturer’s  specifications.  For  this  problem,  acceleration  is  being 

fin'' 


expressed  in 
sown  in  Table  G.2. 


^  so  these  constants  are  converted  to  those  units.  These  are 


Table  G.2.  Final  Accelerometer  Calibration  Constants 


Accelerometer  Serial  ^ 
AH  72 
BA  26 


Calibration  Constant  (mV/in/s2) 
2.28 
2.58 


The  calibration  procedure  for  the  force  gauge  used  is  somewhat  more 
complex,  although  the  idea  behind  it  is  very  simple.  The  problem  is  only 
complicated  slightly  because  it  is  necessary  to  test  the  force  gauge  as  it  is  to 
be  used  on  the  impact  hammer.  There  are  two  major  reasons  why  this  must 
be  so: 

1 )  The  impact  hammer  has  a  sd't  tip  attached  in  order  to  better  excite 
the  low  frequency  modes  that  are  of  interest.  This  soft  tip  has  nonlinear 
characteristics  that  must  be  accounted  for  that  are  difficult  to  predict 
analytically. 
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2)  The  impact  hammer  behaves  differently  depending  on  the  way 
different  users  strike  the  test  specimen  with  it.  This  is  due  to  factors 
such  as  the  length  and  mass  of  the  user's  arm,  among  other  things. 


To  calibrate  the  force  gauge,  an  accelerometer  is  attached  to  a  small 

object  of  known  mass,  suspended  as  a  pendulum.  The  accelerometer  is 

mounted  on  one  side  of  the  mass  so  that  its  sensitive  axis  is  aligned  with  the 

pendular  motion.  The  mass  is  struck  on  the  other  side  with  the  impact 

z(s) 

hammer,  and  a  transfer  function  of  acceleration  to  force,  H(s)  -  ,  is 

computed.  From  Newton’s  Second  Law  it  known  that  this  quantity  should 
equal  the  reciprocal  of  the  known  mass  (including  the  mass  of  the 
accelerometer).  If  the  same  force  gauge  and  accelerometer  are  to  be  used  to 
perform  the  experimental  modal  analysis,  they  may  be  calibrated 
simultaneously  this  way  with  a  single  calibration  constant  calculated  for  the 
problem.  For  this  test,  ail  measurements  will  be  made  relative  to  the  AH72 
accelerometer,  so  the  force  gauge  will  be  calibrated  with  respect  to  that. 

The  response  varies  slightly  with  frequency,  and  an  average  value  is 
calculated  for  the  frequency  range  of  interest.  The  final  values  are  shown  in 
Table  G.3. 

Table  G.3.  Calibration  Constant  for  Force  Gauge 

H(s)  Calibration  Constant 

(V/in/s2)  (V/lb) 

114.6  0.278 


JL 

m 


(  S2  \ 


ylb-iny 


0.002425 
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Data  Acquisition 

Now  that  the  system  is  calibrated  it  is  possible  to  measure  transfer 
functions  for  several  points  along  the  beam.  After  acquiring  the  data  on  the 
B&K  2032,  it  is  transfered  to  the  StarModal  package  for  further  processing. 
This  package  curve  fits  the  data  and  estimates  the  modal  parameters  of  each 
measurement.  All  measurements  are  taken  with  the  actuator  attached  but 
unpowered.  All  measurements  are  taken  by  striking  the  impact  hammer  at 
the  point  where  the  actuator  is  attached,  13.  The  results  of  this  for 
measurements  taken  at  the  two  points  of  interest  are  shown  in  Table  G.4. 

Table  G.4.  Modal  Data  From  Experimental  Measurements 


Meas.  Pt 

x5 

x5 

xll 

xll 

Mode 

1 

2 

1 

2 

Magnitude 

3.98 

121.9 

18.95 

144.83 

Phase  (Deg) 

19.1 

0.363 

0.597 

190.6 

Damping  (%) 

6.44 

2.04 

6.66 

2.03 

fn  (Hz) 

6.24 

34.7 

6.25 

34.7 

Data  Analysis 

The  final  task  remaining  is  to  extract  the  mode  shapes  from  the  modal 
data  found  above.  The  mode  shapes  are  related  to  modal  data  by  the 
following  relation  (12: 195-198), 
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I  Sr  Ijk  ■  I  Gjk(i  CJf )  j  2*^r 


(G.l) 


where 

£r  is  the  rth  mode  shape  of  the  structure 

%  is  the  natural  frequency  of  the  rth  mode 

^  is  the  damping  coefficient  of  the  rth  mode 

Gjk  is  the  magnitude  of  the  transfer  function  at  measured 

at  point  k,  and  excited  at  point  j  (sign  is  determined  from  phase 

information) 

The  resulting  normalized  mode  shapes  are 


As  a  check  on  the  reasonableness  d*  this  result,  these  numbers  compare 
favorably  with  those  found  from  the  22  DOF  finite-element  model. 
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Gain [a_,b_] 

Blcx:k[ Cshat,rc,rcl,rpl,rp2) , 

{4483.2,-11815,0,0,0,0,0,0}}; ■ 

■>f  em-  {0,0,0,0,0,0,124.6,-93.041}; 

K-{)cl,k2,k3,k4,)c5,k6,k7,k8}; 

wn“{40 . 046, 218 . 73 }  ; 

poles-  {(-a+Sqrt[l-a  a]  I)*vm[[l]], 

(-a-Sqrtil-a  a]  I)*wn[[l]l, 

(-b+Sqrt[l-b  b]  I)*wn{{2]], 

(-b-Sqrt[l-b  b]  r)*wn{{2]]}; 

Sya-(teta4Lt{  [1]  ] )  (beta-shat[  [2]  ] )  (beta-shat{  [3]  ] )  ( beta-shat [  [4]  ] )  ; 

alpha-Det  (beta  identityMatrix{8}-af6m+0uter[Tajnes,bfem,k]] , 
alpha-Collect[alpha,beta] ; 
omega-Collect [ omega , beta ] ; 

quot-PolynomialQuotient {alpha, omega, beta}  ; 

quot-Collect ( quot , beta ] ; 

rem-PolynomialRemainder (alpha , omega, beta] ; 

rem-Collect  { rem ,  beta  ]  ,• 

rc-CoefficientList[quot,beta] ; 

rpl-rc[[4]]  rc[[3]]-rc[[5]]  rc[[2]]; 
rp2-rpl  (rc((2]]-rc[{4]]  rc[(l]l); 
rcl-Append[rc,rpl] ; 
rc2-Append [ rc 1 , rp2 ] ; 
remout^CoefficientULst  [reiH/beta]  / 

equal -Solve ( {remout-- { 0 , 0 , 0 , 0 }}, {k3 , k4 , k7 , k8} J ; 

equality-Simplify [equal] ; 

equal2-Solve( {remout““{0,0,0, 0} } , Ck5,k6,k7,k8}] ; 

S^-{{0^0?i?0}!’{o!o?0?i}!{-37268,8235.8,0,0},{44832,-11815,0,0}}; 

)red- {0,0, 121. 61, -79. 026}; 

kred-{Kl,K2,K3,K4};  '  ^ 

roots-Det[s  iden tityMatrix [4 ] -ared+Outer[ Tunes, bred, kred] ]  ; 
sroots-Simplify[roots]; 

listroots-CoefficientList[sroots,s] ; 

Expoles-(s-poles{[l]] )  (s-poles{ [2] ] )  {s-poles[ [3] ] )  (s-polesI[4]]); 

ExpPoles-CoefficientList [Expoles, s] 

Solution-Solve[  {listroots"ExpPoles} ,  {K1,K2,K3,K4}]  ] 
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